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Summary 

Extensive slow-crack-growth (SCG) analysis was made using a primary exponential crack- velocity 
formulation under three widely used load configurations: constant stress rate, constant stress, and cyclic 
stress. Although the use of the exponential formulation in determining SCG parameters of a material 
requires somewhat inconvenk nt numerical procedures, the resulting solutions presented gave almost 
the same degree of simplicity in both data analysis and experiments as did the power-law formulation. 
However, the fact that the ineit strength of a material should be known in advance to determine the 
corresponding SCG parameters was a major drawback of the exponential formulation as compared w ith 
the power-law formulation. 


Introduction 

Advanced ceramics are candidate materials for structural applications in advanced heat engines and 
heat recovery systems. The major limitation of these materials in hostile environments, particularly at 
elevated temperatures, is slow -crack-growth (SCG)-associated failure, where slow' crack growth of 
inherent defects or flaws can occur until a critical size for catastrophic failure is reached. To ensure 
accurate life prediction of ceramic components, it is important to accurately evaluate the SCG parameters 
of a material with specified loading and environmental conditions. 

Life prediction (or SCG) parameters of a material depend on what type of crack-velocity formulation 
is used to determine them. The power-law crack-velocity formulation has been used for several decades 
to describe SCG behavior of a variety of brittle materials ranging from glass and glass ceramics to 
advanced structural ceramics. The main advantage of the power-law formulation over other crack- 
velocity formulations lies in the simplicity in its mathematical expression for lifetime analysis. It has 
also been observed that the power-law formulation has described adequately the SCG behavior of many 
brittle materials. Because of these merits, the power-law formulation has been used in two recent ASTM 
test standards (refs. 1 and 2) to determine SCG parameters of advanced ceramics in constant stress rate 
testing at both ambient and elevated temperatures. Alternative crack-velocity formulations take exponen- 
tial forms to account for the influence of other phenomena (such as a corrosion reaction, diffusion control, 
thermal activation, etc.). However, these exponential forms generally do not result in simple mathematical 
expressions of life prediction formulation, although the forms might better represent the actual SCG 
behavior of some materials. Because of this mathematical inconvenience, the exponential crack-velocity 
formulation has rarely been u ed for brittle materials as a means of life prediction methodology in testing 
or analysis. 
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In this report, the exponential crack-velocity formulation was analyzed to achieve a more convenient 
and simplified life prediction analysis compared with the previous exponential crack-velocity-based 
analyses. The numerical analysis presented here was made for three widely utilized load configurations: 
constant stress rate (dynamic fatigue), constant stress (static fatigue or stress rupture), and cyclic stress 
(cyclic fatigue). The resulting analysis obtained with the exponential formulation was compared with that 
of the power-law formulation to assess which would yield a better life prediction methodology in terms 
of accuracy and convenience in testing and analysis. To the authors’ best knowledge, no analytical study 
on slow crack growth has been done previously using the exponential formulation under cyclic loading. 

In the following reports (parts 2 and 3 of this series) the merits and limitations of the exponential formula- 
tion will be further described in detail using a variety of SCG data determined for many glasses and 
advanced ceramics at both ambient and elevated temperatures. 

All symbols used in this report are listed in the appendix. 

This work was sponsored in part by the HOT/PC and the ZCET projects at the NASA Glenn Research 
Center, Cleveland, Ohio. 


Theoretical Background 

Power-Law Formulation 

The widely utilized empirical power-law crack-velocity term for above the fatigue limit is expressed 
in the form (ref. 3) 


da_ 
d t 



( 1 ) 


where 

v crack velocity 

a crack size 

t time 

Kj mode I stress intensity factor 

K JC mode I critical stress intensity factor (or fracture toughness) 

A, n material- and environment-dependent SCG parameters 


Typically, SCG testing to determine related SCG parameters is performed by applying constant stress 
rate, constant stress, or cyclic stress loading to ground-test specimens. Constant stress rate testing deter- 
mines strength as a function of applied stress rate, whereas constant stress and cyclic stress testing mea- 
sure time to failure as a function of applied stress. The strength in constant stress rate and the time to 
failure in constant stress and cyclic stress tests can be analytically derived to give the following relations 
(refs. 4 and 5): 


Of = D d 6 l/{ ' ,+l) (2) 

t fs = D s o~ n (3) 


i 
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t fc ~ D c G max 


(4) 


where Of is the fracture stress corresponding to the applied stress rate o in constant stress rate testing, 
tf s is the time to failure subjected to a constant applied stress o in constant stress testing, and t^ c is the 
time to failure subjected to cy .lie loading with a maximum stress (T max in cyclic stiess testing. The pa- 
rameters represented by D ' s tire expressed as follows (refs. 4 and 5): 


D,,= 


B(n + 1)S''~ 2 


1 , '(;;+! ) 


(5) 


D s = BS" 


( 6 ) 


D c = 



( 7 ) 


where B = 2 K, c l[AY 2 (n - 2)] where Y is the crack geometry factor in the relation K, = Yo a 1/2 ; S j is the 
inert strength at which no slow crack growth occurs; the function f(t ) is a periodic function in cyclic 
loading specified in o(t) = O n ax /( f ) in ;i range of 0 1; and x is the period. The SCG parameters n 

and D (and B or A) can be obtained by a linear regression analysis with experimental data in conjunc- 
tion with the corresponding equation, either (2), (3), or (4), depending on the type of loading. Hence, it is 
straightforward to determine SCG parameters n and D by least-squares fitting of the data, which is the 
most advantageous feature of the power-law crack-velocity formulation. This convenience and merit in 
mathematical simplicity in addition to the use of routine test techniques have led for several decades to 
the almost exclusive use of the power-law crack-velocity formulation in life prediction analysis and test- 
ing for many brittle materials over a wide range of temperatures. 


Exponential Formulation 

Fracture-mechanics-based modeling typically offers a framework in which lifing can be made. 
However, long-term life prediction is sensitive to the relation between the slow crack velocity and the 
stress intensity factor, which depends on many factors itself. As a result, several different exponential 
crack-velocity formulations that have been previously proposed are based on these other factors, which 
include the presence of a chemically assisted corrosion reaction (ref. 6), diffusion-controlled stress 
rupture (ref. 7), a thermally activated process (ref. 8), a chemical reaction with constant crack-tip 
configuration (ref. 9), kinetic crack growth (ref. 10), and others (ref. 1 1). The generalized exponential 
crack- velocity forms thus proposed are 



(8) 


v = A 


£j_ 

\ K ic 


exp 




( 9 ) 
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v = A 


( 10 ) 



v = A exp 


K, 


K 


IC 


( 11 ) 


/ 

v = A 


_Kj_ 

K IC 



Kj 


K 


IC ) 


( 12 ) 


where A and n are SCG parameters and are different from those used in the power-law formulation. 
Unlike the power-law crack-velocity formulation, the exponential crack-velocity forms do not yield 
simple, analytical expressions of either the resulting strength as a function of applied stress rate in con- 
stant stress rate testing or of the resulting time to failure as a function of applied stress in constant stress 
testing or maximum applied stress in cyclic stress testing. Several attempts have been made under both 
constant stress rate and constant stress loading to obtain corresponding lifetime expressions through nu- 
merical integration incorporating linear (refs. 12 and 13) or nonlinear (ref. 14) regression analysis. 
However, this approach still involves complexity in regression technique as compared with the simple 
least-squares approach in the power-law formulation. 

Slow-crack-growth analyses of three load configurations of constant stress rate, constant stress, and 
cyclic stress were made in this section to obtain simpler, representative equations through numerical 
solution, which in turn makes the use of regression analysis easier to determine corresponding SCG 
parameters comparable to the case of the power-law formulation. Trantina (ref. 12) used the exponential 
crack-velocity forms (eqs. (8) to ( 10)) to determine the approximate time-to-failure equations under con- 
stant stress rate and constant stress loading and showed that the coefficients of the exponential equations 
(8) to (10) were insignificant except at very low fracture stress. Ritter et al. (ref. 13) used the exponential 
form (eq. (8)) to determine the relation of strength versus stress rate in constant stress rate loading via a 
numerical method for indentation cracks that possess a residual stress field around the indent. For the pur- 
pose of simplicity and generalization, equation (8) was chosen for the present analysis. This equation was 
taken from Wiederhom and Bolz (ref. 9), who modified the original Hillig and Charles (ref. 6) exponen- 
tial formulation. An additional analysis using other crack- velocity forms (eqs. (9) to (11)) was also made 
and the results will be discussed in the section Other Exponential Formulations. 

To minimize the number of parameters to be specified (such as A, a , a, S f , A^ c , and /), it is con- 
venient to use a normalized scheme, as used previously for the power-law velocity formulation (refs. 15 
to 17): 


K = 


Kj 


K 


IC 


a i 


a i 


a = - 


Si 


. * C 

a = — - 


' max 

Si 


( 13 ) 


where 


K* 

stress intensity factor (SIF) 

r 

time 

C* 

crack size 
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a * applied stress 

applied stress rate 

rc* maximum applied stress 

u max r 1 

in cyclic loading, and ci- is the critical crack size in the inert condition or is the initial ciack size. Using 
these variables, the exponential crack-velocity equation (8) can be normalized as follows. 


dC* _nK' 

Si: ^ 

dr 


(14) 


The corresponding normalized SIF K is expressed in load configurations of constant stiess late, con- 
stant stress, and cyclic sinusoidal stress, respectively, as 

j.'* =d*r*(c *)' 2 ( l5 > 


K 



(16) 


K 



(17) 


where R is the stress (or load) ratio, defined as R = o min /O max , in which o min and C max are the mini- 
mum and maximum applied stresses, respectively, applied in cyclic loading; and co is the angular 
velocity. As typical for ceramics, the crack size at instability in either an inert or fatigue environment 
was assumed to be small compared with the body of the specimens or components (i.e., an infinite-body 
assumption). Differential equation (14) was solved numerically using a fourth-order Runge-Kutta method 
for each respective loading configuration. The initial condition was C = 1.0 at T* =0, and the instability 
conditions were K* — 1 .0 and ti AT /d C > 0. In cyclic loading, the frequency was taken as arbitrary values 
of coa -/A > 10 8 , depending on the values of maximum applied stress and n. The effect of frequency on 
the solution is discussed in the section Cyclic Stress Loading. 


Results of Numerical Analysis 

Constant Stress Rate leading 

The results of the numeri cal solution of normalized fracture stress (strength) Cf as a function of 
normalized stress rate 6 are shown in figure 1 for values of n ranging from 5 to 100. As seen in the 
figure, for a given n, strength decreases with decreasing stress rate and represents the susceptibility to 
slow crack growth. The rate of decrease in strength with decreasing stress rate becomes more significant 
with lower n values, analogous to the case for the power-law formulation, indicating that the lower n 
gives rise to the greater SCG susceptibility and vice versa. The strength approaches its corresponding 
inert strength as stress rate increases to a certain value at which no slow crack growth occurs. Likewise, 
the strength converges close io zero as the stress rate approaches In o = 0. A linear relationship between 
a* f and In a* holds for most of n values within the range of o} = 0.2 to 0.9 with correlation 
coefficients of r 2 > 0.9975. 
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Apparent SCG parameter, n‘ 



-10 0 10 20 30 40 50 60 70 80 90 1 00 


Normalized stress rate, In 

Figure 1 . — Numerical results of normalized strength a/ as a function of normalized stress 
rate &* in constant stress rate loading for different values of SCG parameter n. 





0 20 40 60 80 100 120 0 20 40 60 80 100 120 

True SCG parameter, n True SCG parameter, n 

Figure 2. — Relationship between true SCG parameter Figure 3. — Relationship between intercept p and true 

n and apparent SCG parameter n' in constant stress SCG parameter n in constant stress rate loading, 

rate loading. 
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A linear regression anal> sis of o f and In c in the range of <5 f - 0.2 to 0.9 was made to 
determine the slope and intercept of each individual curve for a given n, based on the following relation. 

a} =-lno + p (18) 

■ n' 

where 1 In' and p are the slo| e and intercept, respectively. A comparison of the true n (an input datum) 
and the apparent (calculated) is shown in figure 2, where a linear relationship between /? and n' is 

evident (except for lower n \alues, particularly for n < 10). Hence, the overall relationship between n 
and n' can be approximated as 


n = 0.9775;? + 1.7384 (19) 

with a correlation coefficient r = 0.9995. Since the difference between ri and n was >8 percent for n < 10 
and <3 percent for n > 20. a further approximation of equation ( 19) can be made for n > 20 as follows: 

n' ~ n (20) 

The relationship between the intercept |3 and n is shown in figure 3. The value of p decreases with 
increasing n and becomes in signif icant, approaching zero, when n > 20. The overall relationship 
between p and n in the range of n — 5 to 100 was 


P = 2.666(») 1-7 > 


( 21 ) 


with a correlation coefficient = 0.9973. For a nonnormalized expression, equation (13) is used to 
reduce equation ( 1 8) to 


a t 1 . ■ 

-f- = — lna + x 
5 / /? 


( 22 ) 


where 


X = 



n \ 


a i 

ASj 


+ P 


(23) 


SCG parameters n ' and '/_ in constant stress rate loading can be obtained from the slope and inter- 
cept by a linear regression analysis of ( oy IS j) versus In c . With n‘ thus calculated, n can be evaluated 
from equation (19). The parameter A can be evaluated using equation (23) from calculated X together 
with P (eq. (21)) and known values of a, and 5,. The solution presented in this study (eq. (22)) is much 
simpler compared with the previous solution by Trantina (ref. 12) in which a simple linear regression 
would hardly be applicable because of the complex functional form of the solution, as shown below 
(note that Of is present in botl sides of the equation): 


°f_ 

Si 



(24) 
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Therefore, when determining SCG parameters, the current solution (eq. (22)) significantly eliminates the 
complexity associated with a regression analysis that would be encountered in the previous solution (eq. 
(24)). For the case of n > 20, based on the results of figures 1 and 3 for <j f > 0.4, P is negligible (with a 
maximum of about 7 percent) compared with a f. Hence, equation (21 ) reduces to 


p-o 


(25) 


which results in fin (a i /AS i )]/«. Likewise, in this case, n ' - n from equation (20). 

A distinct difference in functional expression between the power-law and exponential formulations 
is that in the power-law formulation, log Gf is plotted as a function of log d, whereas in the exponential 
formulation, Of/S i is plotted as a function of In c (fig. 1 shows this for the normalized parameters, a f as 
a function of In d ). Hence, the knowledge of inert strength in the exponential formulation is a prerequi- 
site to determining SCG parameters ri and which is a disadvantage compared with the case (eq. (2)) 
of the power-law formulation. Note that the power-law formulation does not require any prior know ledge 
of inert strength to determine SCG parameters n and D. Although not presented here, the numerical re- 
sult was plotted as log a^as a function of log a for different n values in the same way that is used for 
the power-law formulation. The resulting plots, however, showed appreciable nonlinearity, which made 
linear least-squares fitting inapplicable in determining the related SCG parameters. 


Constant Stress Loading 

The numerical results of normalized time to failure Tj as a function of normalized applied stress 
a* are shown in figure 4 for values of n ranging from 5 to 100. The general trend of the solution can 
be summarized in terms of (1) the convergence of In Tf close to zero with g -► 0, (2) the increased SCG 
susceptibility with decreasing n values, and (3) the linearity between In Tf and G in the range of 
G = 0.2 to 0.9. As a consequence, the relationship between normalized time to failure and normalized 
applied stress within the linear region can be written as 

\nT* f = -/?'c* + p (26) 

zfc jij ry 

The linearity between In Tf and a is manifest when the correlation coefficient r~ > 0.995 for each 
curve is considered. Hence, ri and p can be determined with a reasonable accuracy by a linear regres- 
sion analysis based on the results in equation (26). The relationship between ri and n is shown in figure 
5 and has the following relation: 


ri — 0.9827/? + 3.3440 (27) 

with r 2 = 0.9997. The difference between ri and n was >8 percent for n < 30 and <5 percent for n > 40 
so that a further approximation of equation (27) can be made for the case of n > 40 as 

ri ~n (28) 

The function of P with respect to n is depicted in figure 6, where the intercept P decreases with 
increasing n values, resulting in the best-fit relation 
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Normalized applied stress, a* 

Figure 4.— Numerical results of normalized time to failure Tf as function of normalized 
applied stress <r in constant stress loading for different values of SCG parameter n. 



True SCG parameter, n 

Figure 5. — Relationship between true SCG parameter 
n and apparent SCG parameter n' in constant stress 
loading. 



True SCG parameter, n 

Figure 6. — Relationship between intercept p and true 
SCG parameter n in constant stress loading. 
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(3 = -].913 + 4.985e 


(29) 


- 0 . 049 « 

with r 2 = 0.9907. 

For the nonnormalized expression, equation ( 13) is used to reduce equation (26) to 

\nt f =-n'^- + x (30) 

where 

X = ln^j+p (31) 

Hence, ri and % in constant stress loading can be obtained from the slope and intercept, respectively, by 
a simple linear regression analysis of the data In ^ as a function of a/S,-. With ri calculated, n can be 
evaluated from equation (27). The parameter A can be estimated from equation (3 1 ) with calculated % 
together with p (eq. (29)) and known values of a } . The solution presented here is much simpler than the 
rather complex one proposed by Trantina (ref. 12), in which the slope of the relation In t y. versus <5iS { was 
n + (S^/a). Both analyses, however, would give a similar result when n is significantly greater than S 7 /g. 

A notable difference in constant stress loading analysis between the power-law and exponential for- 
mulations is that in the power-law formulation, log /y is plotted as a function of log G as seen in equation 
(3). However, in the exponential formulation, In is plotted as a function of <y/S } . Hence, as in the case 
of constant stress rate loading, inert strength must be known to determine ri and which is a distinctive 
drawback of the exponential formulation as compared with the power-law formulation. Unlike the 
power-law formulation (eq. (3)), there was significant nonlinearity in plots of log Tf as a function of 
log a" (not shown here), which made linear least-squares fitting inapplicable to determine the related 
SCG parameters. 

Cyclic Stress Loading 

The results of the numerical solution of normalized time to failure Tf as a function of normalized 
maximum applied stress a„ iax in cyclic sinusoidal loading with two stress ratios R of 0. 1 and 0.5 are 
shown in figure 7 for values of n ranging from 5 to 80. Similar to the case of constant stress loading, 
this plot is linear with respect to G max in the range 0.2 to 0.9 and converges close to zero with a further 
decrease in a max . The linearity between In Tf and G max was evident considering the correlation coef- 
ficient of r 2 > 0.997. Also, note that the effect of the R- ratio on the solution for a given n value is insig- 
nificant. Based on the results in figure 7, similar to the case of constant stress loading, the relationship 
between In Tf and G max can be described as follows: 

In Tf = -/fa max +p (32) 

where ri is the slope and (3 is the intercept, which can be determined from the numerical results using a 
linear regression analysis based on equation (32). 

Figure 8 shows the relationship between ri and n for R- ratios of 0.1 and 0.5. For comparison, the 
result determined in constant stress loading (i.e., R = 1.0) from figure 5 was also included. Note that the 
SCG analysis in cyclic stress loading reduces to that of constant stress loading when R = 1.0 (also seen in 


NASA/TM — 2002-211 1 53/PART 1 


10 


20 r- 


SCG 

parameter, 
n I 



-120 I 1 1 * 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Maximum normalized applied stress, wmax 

Figure 7.— Numerical results of normalized time to failure Tf as function of maximum 
normalized applied stress a^ ax in cyclic (sinusoidal) stress loading for R-ratios of 0.1 
and 0.5 and diffc rent values of SCG parameter n. 


eq. (17)); hence, constant sirevs loading can be regarded as one of the generalized cyclic loading configu- 
rations. A good linear relatioi between /i’ and n was found for both /? = 0. 1 and 0.5 with r~ > 0.999; 


n ' = 0.9777/r + 2.5296 for R = 0.1 

n' = 0.9772/? + 2.5411 for/? = 0.5 (33) 

As seen in figure 8, no significant difference in the n’-n relation exists for R - ratios ranging from 0.1 to 1.0 
(constant stress loading): this relationship is also observed in the power-law formulation (refs. 5, 17, and 
18). The difference between .<?' and // was >7 percent for // < 20 and <3 percent for n > 40, so that a fur- 
ther approximation of equation (33) can be made for R = 0.1 and 0.5 for the case of n > 40 as follows: 

n' ~ n (34) 

The relationship between the intercept (3 and n is depicted in figure 9, where the result from figure 6 
for constant stress loading </? = 1.0) was also included for comparison. The intercepts for R = 0. 1 and 0.5 
decrease more monotonically w r ith increasing n values than that of R = 1 .0. The best-fit equation, similar 
to equation (29) in constant stress loading, was obtained for each /?-ratio: 

p = 0. 1 409 + 3.559e -() 0737 " for R = 0. 1 

(^O.l^-^^e -0 ' 0857 " for/? = 0.5 (35) 

with the correlation coefficient of r 2 > 0.991 . 
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Stress 

ratio, 



True SCG parameter, n 

Figure 8. — Relationship between true SCG parameter 
n and apparent SCG parameter n' in cyclic (sinu- 
soidal) stress loading with R-ratios of 0.1, 0.5, and 
1 . 0 . 


Stress 

ratio, 

R 



0 20 40 60 80 100 120 


True SCG parameter, n 

Figure 9. — Relationship between intercept p and true 
SCG parameter n in cyclic (sinusoidal) stress load- 
ing with R-ratios of 0.1 , 0.5, and 1 .0. 


For the nonnormalized expression, equation (13) can be used to reduce equation (32) to 


lnt f =-n'^ + X 


where 



(36) 


(37) 


Therefore, n f and x °f the exponential formulation in cyclic stress loading for a given 7?-ratio can be 
obtained from the slope and intercept by a linear regression analysis of the data of In y- as a function of 
G max /Sj* With n’ calculated,/* can be evaluated from equation (33). The SCG parameter A can be 
estimated from equation (37) with calculated P (eq. (35)), and known values of a i 

In cyclic stress loading, a distinct difference between the power-law and exponential formulations is 
that in the power-law formulation, log yds plotted as a function of log a max as seen in equation (4) 
whereas in the exponential formulation, In *yis plotted as a function of C max /S i . Hence, as in the cases of 
constant stress rate and constant stress loading, the inert strength of a material must be known in advance 
for cyclic loading to determine the corresponding SCG parameters, a clear disadvantage of the exponen- 
tial formulation compared with the power-law formulation. Although not presented here, it was found that 
in plots of log T j as a function of log o max , typical of the power-law formulation (eq. (4)), there was a 
considerable nonlinearity that made linear least-squares fitting inapplicable in determining the related 
SCG parameters. 

Typical examples of the effect of frequency on the time to failure in cyclic stress loading are shown 

* 

in figure 10 for /?-ratios of 0. 1 and 0.5 with n = 20 and o max = 0.9. The number of cycles to failure Nj is 
calculated using the relation 
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Fig jre 10. — Normalized time to failure Tf as function 
cr number of cycles to failure Nf in cyclic (sinu- 
soidal) stress loading for R-ratios of 0.1 and 0.5 
under for n = 20 and o^ax = 0.9. Arrow indicates 
n jmber of cycles corresponding to unity. 


where (0 a-JA is the input value for numerical procedures, as shown in equation (17). As seen from 
figure 10, Tf decreases rapid y around N f = 10“' to 10° and tliereafter approaches a plateau that corre- 
sponds to the exact solution oi time to failure. In other words, as long as the total number of cycles to 
failure is > 1 , the solution of ti me to failure is converged and is thus independent of either the number 
of cycles or frequency. Hence in the numerical procedure, the value of coa-M was chosen to fulfill this 
requirement and started with a minimum value of 10 8 , depending on Tf . This frequency independency 
in the exponential formulation is the same as that in the power-law formulation (refs. 5 and 17). For the 
power-law formulation, the time to failure remains unchanged with frequency, which is attributed to the 

T 

functional form of (1 / x)J [/(r )]"dr in equation (7) (ref. 5). It can be easily shown that the period t of 

o 

any applied loading cycle is a ways cancelled out in the formulation. It should be noted that in this analy- 
sis, the presence of another damage mechanism (refs. 19 to 21 ) in addition to slow crack growth was not 
assumed to occur in cyclic loading. Only slow' crack growth was considered as the unique mechanism 
leading to the failure of a material. 


Comparison of Constant Stress and Cyclic Stress Loading Lifetimes: the //-Ratio 

The ratio of constant stress and cyclic stress loading lifetimes, the //-ratio, with a condition of c in 
constant stress loading equal to a max in cyclic stress loading (a* = a^ ax ), has been frequently used in 
the power-law' formulation (ri fs. 5, 17, and 18) to quickly compare lifetimes of constant stress and cyclic 
stress loading. As done customarily for the power-law formulation, the //-ratio was also calculated for the 
exponential formulation. The //-ratio is defined as (refs. 5, 17, and 18) 

h (39) 


NASA/TM — 2002-211 153/PAR 11 


13 



Stress 

ratio, 



Figure 1 1 . — Ratio of constant stress to cyclic stress 
loading lifetimes, h-r atio, as function of SCG param- 
eter n for different fl-ratios in cyclic (sinusoidal, 

°max = 0.7) loading. Each line represents best fit 

where again tj s and t^, are times to failure in constant stress and cyclic stress loading, respectively. 
Based on the numerical results of time to failure for constant stress and cyclic stress loading, the /7-ratio 
was calculated using equation (39), and the results are presented in figure 1 1 as a function of n for differ- 
ent values of the /?- ratio. The /?-ratio decreases with increasing n, and the rate of decrease with increas- 
ing n is almost the same regardless of /f-ratio up to 0.9. Also, for a given 77, the / 7 -ratio increases with 
increasing /?-ratio. Note that the R- ratio of 1.0 represents the case for constant stress loading so that when 
R = 1.0, the numerical solution in cyclic stress loading should reduce to the case of constant stress load- 
ing. This presents another way to check the accuracy of the cyclic stress loading analysis. The / 7 -ratio var- 

+ * 

ies slightly by a factor of 2 between a max = 0.2 and 0.9: the lower o max value gives the higher /?-ratio 
and vice versa. Hence, for a conservative estimate, the higher value of a max = 0.7 was used for the calcu- 
lation of the / 7 -ratio. The maximum difference in life between cyclic stress and constant stress loading is 
approximately 1 order of magnitude, which occurs for R = 0. 1 and n > 80. A similar trend in the / 7 -ratio 
can also be observed in the power-law formulation, as shown in figure 12. Here the plots of the / 7 -ratio 
as a function of n are shown for a range of R- 0.0 to 1 .0, calculated previously for the power-law 
formulation (ref. 17). Unlike the exponential formulation, no effect of a max on the /i-ratio for a given 
/?- ratio had been observed for the power-law formulation. Although no significant difference exists in 
the plots of / 7 -ratio versus n between the exponential (fig. 1 1) and power-law (fig. 12) formulations, 
the overall magnitude of the / 7 -ratio is about 20 percent greater in the exponential than in the power-law' 
formulation. 


Other Exponential Formulations 

A comparison of solutions from other exponential SCG formulations under three loading configura- 
tions is shown in figure 13. The figure presents the results of three exponential formulations of equations 
(9) to (1 1) for n = 20 and 40 and compares them with those of the primary formulation of equation (8). 
The difference in solution between this primary equation and two other equations ((9) and (10)) was 
insignificant, particularly at higher stress rates (in constant stress rate loading, fig. 13(a)) and higher 
applied stresses (in constant stress and cyclic stress loading, figs. 13(b) and (c)), giving rise to a reason- 
able linearity between the dependent/independent variables related. This insignificant difference in solu- 
tion as well as the linearity allows one to conclude that the primary equation would be representative of 
all three exponential SCG formulations considered. By contrast, the remaining second-order formulation 
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Normalized strength, <rf 


Stress 



Figure 12. — Ratio of constant stress to cyclic stress loading 
lifetimes, h-ratio, as function of SCG parameter n for differ- 
ent ft-ratios in cyclic (sinusoidal) stress loading, calculated 
with power-law formulation (ref. 17). Each line represents 
be; t fit. 




Normalized applied stress, o-* 



.120 ^ 1 1 1 1 1 

0.2 0.4 0.6 0.8 1.0 1.2 


Maximum normalized applied stress, cr^ax 

Figure 13. — Results of numefical solutions using three exponential formulations of equations (9) to (11) com- 
pared with the primary exponential formulation of equation (8) for selected SCG parameters of n = 20 and 
80. (a) Constant stress-rate loading, (b) Constant stress loading, (c) Cyclic stress loading with stress ratio 
R = 0.1 . 
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(eq. (1 1)) showed an appreciable deviation and a notable nonlinearity. Therefore, the determination of 
corresponding SCG parameters in this case differs from that of the primary equation and should only 
be attempted under appropriate circumstances so that a simple linear regression could be applied with 
reasonable accuracy. It is expected that considering its functional form, equation (12) also would yield 
results similar to those of equation (11). 

Although the use of the exponential formulations to determine the SCG parameters of a material 
requires somewhat inconvenient numerical procedures, the resulting solutions given in this report would 
have almost the same degree of simplicity in both data analysis and experiments as the power-law 
formulation in constant stress rate, constant stress, or cyclic stress loading configurations. However, 
the knowledge of inert strength of a material should be known beforehand so that the corresponding SCG 
parameters (particularly /?) can be determined, which could be a major drawback of the exponential for- 
mulation. In parts 2 and 3 of this report series, a variety of experimental data from various glasses and 
advanced ceramics at both ambient and elevated temperatures will be used to verify the solutions given 
herein. 


Conclusions 

Based on the numerical solutions of life prediction parameters obtained with exponential formula- 
tions, the following conclusions were made: 

1. In constant stress rate (dynamic fatigue) loading, slow-crack-growth (SCG) parameters can be 
determined by a linear regression analysis of the data of (fracture stress/inert strength) as a function of 
applied stress rate, oy/5- versus In <T, together with the appropriate relations provided. 

2. In constant stress (static fatigue or stress rupture) and cyclic stress (cyclic fatigue) loading, the 
corresponding SCG parameters can be evaluated by a linear regression analysis of the data of time to fail- 
ure as a function of (maximum applied stress/inert strength), In tj- versus c max /S;, in conjunction with the 
pertinent relations provided. 

3. No frequency effect on life and no dependency of SCG parameter n on the /?-ratio, the ratio of 
minimum to maximum applied stress, was observed in cyclic stress loading, much the same as that 
observed for the power-law formulation. The difference in the ratio of constant stress to cyclic stress load- 
ing lifetimes, the //-ratio, between the exponential and power-law formulations was minimal with a maxi- 
mum difference of about 20 percent. 

4. While the numerical solutions using the exponential formulation require somewhat inconvenient 
numerical procedures, they provide almost the same level of simplicity in both data analysis and experi- 
ments as the power-law formulation. However, requiring the knowledge of the inert strength of a material 
to determine corresponding SCG parameters (particularly n) would make the exponential formulation 
more difficult to use in comparison with the power-law formulation. 

5. There is no appreciable difference in solutions between the primary exponential SCG equation 
(used in this analysis) and two other exponential expressions analyzed in this study, so the primary equa- 
tion would be considered representative of all the first-order exponential formulations. 
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Appendix — Sy m bols 


A slow-crack-growth parameter defined in equations ( 1 i and (8) 
a crack size 

B slow-crack-growth parameter, B = 2A^/ [AY 2 (?i - 2)| 

C crack size in normalized scheme of references 15 to 17 

D slow-crack-growth pi rameter defined in equations (5 ) to (7) 
f(t) periodic function, cy< lie loading 

h ratio of constant to cvelic stress loading lifetimes 

K stress intensity factor 

N number of cycles 

n slow-crack-growth p; rameter defined in equations ( 1 i and (8) 

R stress ratio 

n correlation coefficient 

S strength 

T time in normalized scheme of references 15 to 17 
t time 

v’ crack velocity 

Y crack geometry factor 

p intercept of curve in linear regression analysis defined in equations (18), (26), and (32) 

X slow-crack-growth p; rameter defined in equations (23), (31), and (37) 
a applied stress 

applied stress rate 
i period 

(0 angular velocity 
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Subscripts: 

C critical 

c cyclic stress 

d constant stress rate 

/ fracture 

/ mode I 

i inert or initial condition 

max maximum 

min minimum 

5 constant stress 

Superscripts: 

* normalized 

apparent (calculated) 
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